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ABSTRACT 

An improved formula for the timing of binary pulsars that accounts for the relativistic 
deflection of light in the gravitational field of the pulsar's companion is presented, and 
the measurability of this effect together with its variance estimates are discussed. The 
deflection of the pulsar's beam trajectory in the gravitational field of its companion 
leads to variation in the pulsar's rotational phase. This variation appears as the narrow 
sharp growth of the magnitude of the post-fit pulsar timing residuals in the vicinity of 
the moment of the superior conjunction of the pulsar with its companion. In contrast 
to the relativistic Shapiro effect the amplitude of the effect of gravitational deflection 
of the pulsar radio beam has two peaks with opposite signs, which become sharper as 
the inclination i of the pulsar's orbit approaches 90°. The effect under consideration 
influences the estimation of parameters of the relativistic Shapiro effect in the binary 
pulsars with nearly edgewise orbits. Its inclusion in the fitting procedure provide more 
careful measurement of sine of the orbital inclination i as well as the masses of the 
pulsar m p and its companion m c . This permits an improved testing of alternative 
theories of gravity in the strong field regime. Moreover, the measurement of this effect 
permits independent geometric constraints on the position of the pulsar rotational 
axis in space. The effect of the gravitational deflection of light has been numerically 
investigated for binary pulsars with nearly edgewise orbits. It is shown that the effect 
is observed in general only when cosi < 0.003. This estimate becomes less restrictive 
as the pulsar's spin axis approaches the line of the sight. 

Key words: gravitation - pulsars - relativity - stars: neutron - stars: individual (PSRs 
B1855+09, B1913+16, B1534+12) 



1 INTRODUCTION 

High-precision timing observations of binary pulsars on rel- 
ativistic orbits provide unique opportunities for testing of 
General Relativity (GR) through two mechanisms: the emis- 
sion of gravitational radiation and the presence of strong- 
field effects. The effect of emission of gravitational waves 
have been confirmed using the PSR B1913+16 long term 
data with an accuracy better than 0.5% (Taylor & Weis- 
berg 1989; Damour &Taylor 1991). This was accomplished 
by the measurements of three relativistic effects: 1) the 
rate of apsidal motion w, 2) the Doppler and gravitational 
red shifts 7 of pulsar rotational frequency, and 3) the sys- 
tem orbital period decay F%, which is due to the effect of 
the gravitational waves emission, which has been rigorously 
justified from the mathematical point of view by Damour 
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(1983a,1983b) and, independently, by Grishchuk & Kopeikin 
(1983,1986), as well as Schafer (1985). However, timing of 
PSR B1913+16 is not sufficiently sensitive to the strong 
gravitational field effects. And, indeed, the theoretical inves- 
tigations of Damour and Esposito-Farese (1992) have shown 
that some alternative to the GR theories of gravity, which 
include the additional scalar fields, also fit the above men- 
tioned u) — 7 — Pb test, and, hence, can not be distinguished 
from the GR theory. Realizing this fact has led to the devel- 
opment of the parameterized post-Keplerian (PPK) formal- 
ism, which was designed to extract the maximum of possible 
information from the high-precise observational data for bi- 
nary pulsars (Damour & Taylor 1992). The application of 
this methodology to the four pulsars B1913+16, B1855+09, 
B1534+12 and B1953+29 has already ruled out a number 
of the alternative gravity theories (Taylor et al. 1992). 

The PPK formalism is especially effective in discrim- 
inating alternative theories of gravity in the strong field 
regime when there is a possibility of independently measur- 
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ing two post-Keplerian parameters chaacterizing the range 
r = Grric/c 3 and the shape s = ami of the famous Shapiro 
delay in propagation of radio signals of pulsar in the gravi- 
tational field of its companion. The most favorable case for 
the measurements of r and s is when the pulsar's orbit is 
inclined to the plane of the sky at the angle i close to 90°. 
Timing observations of the " non-relativistic" binary pulsar 
B1855+09 by Ryba & Taylor (1991), and by Kaspi et al. 
(1993) have provided the very first such independent mea- 
surements of the parameters r and s that allowed the deter- 
mination of the neutron star mass with a precision better 
than about 18%. It is worthwhile to stress that such mass 
determination does not require a knowledge of the rate of 
relativistic apsidal motion Co and the red shift parameter 7 
as it takes place in the case of the PSR B1913+16 system. 

At the same time when a binary pulsar orbit is viewed 
nearly edge-on (as for B1855+09 and PSR B1534+12) the 
complementary effect in the propagation of the pulsar's ra- 
dio signal may be important. This is due to the well known 
effect of the deflection of light rays in the gravitational field 
of a self-gravitating body (in the gravitational field of the 
pulsar's companion in the event of the binary pulsar system). 
Due to the gravitational deviation of the pulsar's beam from 
a straight line, the observer will see the pulsar's pulse only 
when the pulsar is rotated by an angle compensating this 
deviation. In contrast to the relativistic Shapiro effect, the 
effect of gravitational deflection of light directly contributes 
to the rotational phase of the pulsar (as, for instance, the 
effect of aberration of pulsar's beam does), but this is so 
small that it can be misinterpreted as an additional delay 
in the time-of-arrivals (TOA) of the pulsar pulses. Similar 
to the Shapiro delay, the deflection of the pulsar's beam by 
the gravitational field of its companion manifests itself in 
the timing as the rapid, sharp growth of the magnitude of 
the post-fit residuals of TOA on a short time interval in the 
vicinity to the moment of superior conjunction of the pulsar 
and its companion. The effect in question is superimposed 
on the Shapiro effect, and, in principle, should be explicitly 
taken into account to avoid an incorrect determination of 
the parameters r and s. 

Another important feature of the effect of gravitational 
deflection of light in a binary pulsar is a strong dependence 
of the shape of TOA residuals on the spatial orientation of 
the pulsar's rotation axis (see the eqs. (7), (9), (10) in the 
text below). Thus, the measurement of this effect would help 
to get further constrains on the position of the pulsar's spin 
independently of the polarization observations. Let us em- 
phasize that the effect of the relativistic aberration caused 
by the orbital motion of pulsar in binary systems (Smarr 
& Blandford 1976) also strongly depends on the spin's ori- 
entation, but can not be observed directly due to the same 
time- dependence on the pulsar's orbital phase as the classi- 
cal Roemer effect (Damour & Deruelle 1986). 

Epstein (1977) was, probably, the first who called at- 
tention to the phenomenon of deflection of the pulsar ra- 
dio beam in the gravitational field of its companion when 
cosi < 10 -3 . However, he did not make any calculations 
confirming this estimate. Narayan et al. (1991) also men- 
tioned this effect as a useful tool for investigation of black 
hole - neutron star binary systems and verification of the 
existence of the black hole to an unparalleled degree of cer- 
tainty. Schneider (1989, 1990) calculated the corresponding 



gravitational deflection time delay and lensing amplification 
factor of the pulse intensity in the edge-on binary pulsar, 
but, unfortunately, his presentation is not easily followed. 

In the present paper we try to solve the problem more 
simply, giving both a complete theoretical treatment and a 
discussion of the measurability of the effect of the gravita- 
tional deflection of radio beam in binary pulsars which we 
call briefly the bending delay (BD) effect. It is shown how 
to include this effect in the fitting procedure of TOA for the 
evaluation of the pulsar parameters of binary pulsars whose 
orbital inclinations to the plane of the sky are extremely 
close to 90° . 

In Sec. 2 of the paper we present an improved formula 
for the timing of binary pulsars with the BD effect taken into 
account. In Sec. 3 we investigate, using numerical simula- 
tions, the measurability of the masses and orbital inclination 
in binary systems when the BD effect is taken into account. 
A mathematical procedure for the direct measurement of 
parameters of the BD effect is outlined in Sec. 4 and applied 
to process the observational data for PSR B1855+09. In Sec. 
5 we present a summary of our results and conclusions. 



2 TIMING FORMULA 

The self-consistent mathematical derivation of a timing for- 
mula for the calculation of pulsar rotational phase should 
be based upon the relativistic theory of astronomical refer- 
ence systems together with the matching procedure devel- 
oped previously by Kopeikin (1988), Brumberg & Kopeikin 
(1989a,b) and Brumberg (1991). The exact timing formula 
must include the kinematic effects of the galactic differen- 
tial rotation, orbital and precessional motions of the pulsar, 
relativistic effect of aberration, quadratic Doppler and red 
shift effects, the Shapiro delay, as well as the deflection of 
the pulsar beam in the gravitational field of its compan- 
ion. In the present paper we ignore the effects of differential 
galactic acceleration and precessional motion of the pulsar 
rotational axis suggesting for simplicity that vectors of the 
angular velocity H and spin S of the pulsar are coincident. 
The timing formula obtained is a slightly modified version 
of that elaborated by Damour & Deruelle (1986). Some ad- 
ditional theoretical aspects of the derivation of the timing 
formula are discussed by Klioner & Kopeikin (1994) and 
Kopeikin (1994, 1995). 

Let r\ and A be the polar coordinates of the pulsar's 
spin axis S, where r\ is the longitude of the pulsar's spin 
axis in the plane of the sky measured from the ascending 
node of the pulsar's orbital plane (0° < r\ < 360°), and A 
is the angle between the pulsar's spin axis and the line of 
sight pointing from the terrestrial observer toward the pulsar 
(0° < A < 180°). The time evolution of the pulsar rotational 
phase 4>{T) is given by the equation: 

0(T)/2tt = N(T) = N + u p T + i*> p T 2 + i^T 3 , (1) 

where No is the initial phase; v v , v v , v v denote the pulsar spin 
frequency and its time derivatives shifted by the Doppler 
factor, and T is the pulsar proper time related to the solar- 
system barycentric arrival time t according to the relation- 
ship: 

t - to = T+ A R {T) + A E (T) + As (T) + A A (T) + A B (T) . (2) 
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Here to is the initial barycentric epoch of observation; A R , 
Ae and A a correspond to the well-known Roemer, Ein- 
stein and aberration time delays, whose exact mathematical 
expressions can be found in the paper of Damour & Tay- 
lor (1992); As is the Shapiro delay of radio pulses in the 
gravitational field of companion of the pulsar (Epstein 1977; 
Damour & Taylor 1992): 

A s = -2rlnL, (3) 

L — 1 — ecosit — s{sintj(cosii — e) + (l — e 2 ) 1//2 cos u sin u}, (4) 

which depends on the just introduced post-Keplerian pa- 
rameters r and s as well as the classical orbital parameters: 
eccentric anomaly u, eccentricity e, and longitude of the pe- 
riastron lu. 

Ab in the equation (Q) is the additional bending delay 
caused by the deflection of the pulsar beam in the gravi- 
tational field of its companion. It has been calculated from 
the condition of emission of the pulsar's pulse toward the 
observer (Kopeikin 1992). This condition is expressed in the 
pulsar's proper reference frame by the equation: 



(N ■ (fiA6)) = 0, 



(5) 



which defines a relation between the number of pulse, which 
has been emitted, and the moment of the emission reckoned 
in the pulsar proper time scale T. Here N denotes the spatial 
components of isotropic unit vector tangent to the trajectory 
of the emitted pulse, and directed toward the terrestrial ob- 
server; b is the unit vector in the fiduciary direction of maxi- 
mal pulsar radio emission (it is not necessarily the direction 
of magnetic dipole); Q is the vector of the instantaneous an- 
gular velocity of the pulsar; the signs "•" and "A" denote 
the Euclidean scalar and cross products. We have assumed 
implicitly that the model of the conical rotating beacon for 
the description of the pulsar's emission (Lyne & Manchester 
1988) is true. However, we would like to point out that cal- 
culations being described in this paragraph are still valid in 
the case of a deviation of the polar beam geometry from the 
conical one as it is, for instance, in the model of Narayan & 
Vivekanand (1983). 

Vector N does not keep a constant spatial orientation 
from the point of view of an external observer at rest but 
has small (caused by both classical and relativistic effects) 
spatial variations with respect to the constant unit vector ft 
referred to the rest frame of the binary pulsar bary center. 
To show this we use the simplified version of the relativistic 
relationship between the spatial components of the vectors 
N and ft derived by Klioner & Kopeikin (1992) in the case 
of the relativistic N-body problem: 

N = ft — — (ft A (f p A ft)) + -(n A (v p A ft)) 



2Gm c (ft A (f R A ft)) 



(6) 



c 2 r R r R + (ft - f R ) 

to replace in the equation Q the time dependent compo- 
nents N by the ones of the vector ft. The complete relation 
between N and ft (Klioner & Kopeikin 1992) can be easily 
used for investigations of the second order aberration effects 
and physical parameters of a rotating black hole being con- 
sidered as a possible companion of a pulsar (Narayan et al. 
1991). 

In formula (n), c is the speed of light, r p and f R are 



the radius- vectors pointing respectively from the barycenter 
of the binary system and the center-of-mass of the pulsar's 
companion to the pulsar (do not confuse these vectors with 
the notation r for the parameter characterizing amplitude 
of the Shapiro delay), v p = df p /dT is the orbital velocity of 
the pulsar, and R is the distance between the binary pul- 
sar and the solar system barycenters. Let us note that it is 
possible to equate the vector —ft to the unit vector K (that 
is ft — —K) pointing from the Earth toward the pulsar sys- 
tem as it is usually done, though such procedure requires an 
appropriate justification from a theoretical point of view. 

After substituting expression (6) into equation (5), ex- 
pressing vectors N, Q, and ft through the Euler variables, 
and solving the obtained trigonometric equation with re- 
spect to the fast Euler variable as described in (Doroshenko 
et al. 1995 we get the counting formula (1) with additional 
phase corrections due to the classical and relativistic terms 
in the relationship (6). To be specific, the second term in the 
right-hand side of this formula gives an additional but neg- 
ligibly small classical parallax delay in the rotational phase 
of the pulsar (see, however, the paper of Kopeikin (1995) 
where the influence of the orbital parallax effect on the 
propagation time of pulsar radio signals is shown to be de- 
tectable); the third term in the right-hand side of (6) leads 
to the aberration delay Aa(T) depending on two directly 
unobservable parameters A and B introduced by Damour & 
Deruelle (1986), and the very last term defines after straight- 
forward calculations the gravitational bending delay being 
under consideration: 



A B = L 1 (C sin(w + A e ) + D cos(w + A e )) . 



(7) 



Here the true orbital anomaly A e is related to u by the well- 
known equation: 



A e = 2 arctan 



•l + eV/2 u 
tan- 



(8) 



a R designates the semi-major axis of the relative orbit of 
the pulsar with respect to its companion, and the constants 
C and D have been introduced as the new post-Keplerian 
parameters: 

cr(l-s 2 ^ 2 



c 



D = 



nv p a R 



cost; 
sin A ' 



sin X] 



-Kv v a R sin A 



(9) 



(10) 



where a R — a p (m p + m c )/m c , and a p is the major semi- 
axis of the pulsar orbit with respect to the barycenter of 
the binary system. It is interesting to note that the final 
expression for the BD effect, formula (Q), does not include 
any dependence on the spatial orientation of vector b. The 
same is true for the aberration delay A a- 

For a better understanding of influence of the BD effect 
on the TOA of the pulsar pulses let us compare the behavior 
of the functions As and As by inserting parameters of the 
system PSR B1855+09 into the expressions and (R). The 
explicit dependence on time of the Shapiro (eq. (H)) and 
BD (eq. (0)) effects are shown in Fig.l , which has been 
drawn using observational parameters of the PSR B1855+09 
reported by Kaspi et al. (1994). 

For illustrative purposes we have chosen the artificially 
small numerical value for A = 0°.005, and two different val- 
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arrival times tu with the help of the equations (|l|) - M) and 
(@) " (fiot) - Assuming that the differences 5ej = — e' ' 
between the (yet unknown) true and initial numerical values 
of the parameters are small, one can represent an expression 
for the total derivative of the phase residuals (the differential 
timing formula) in the form: 



R(t) = ^2b,- Se t 



(11) 



Figure 1. The illustrative example of the TOA residuals for 
PSR B1855+09 caused by the Shapiro delay and deflection of the 
pulsar beam in the gravitational field of its companion - smooth 
or dotted curves. For more explanations see the text. 



ues for the angle: r\ = —85° (solid curve, D > 0, C > 0), 
and r\ — 85° (dotted curve, D < 0, C > 0) to depict the 
BD effect in the most clear form. One can see from Fig.l 
that the modulation of TOA due to the BD effect is com- 
pletely different and quite sharper than that caused by the 
Shapiro delay. The total observational timing effect will be 
the sum of the Shapiro and BD delay curves. It is important 
to note that the asymmetry in the curve describing the BD 
effect does not depend on the value of the angle A (since 
<sinA < 1, when 0° < A < 180°) and is sensitive to the 
numerical value of the angle r\ only. Therefore, we can say 
nothing from the investigation of the BD effect whether the 
pulsar spin axis S is along the direction to observer or in 
the opposite one. 

The considered example shows that there is in principal 
a possibility to detect the BD effect in the binary pulsars 
with nearly edgewise orbits. However, this effect was never 
included in existing pulsar timing algorithms. Therefore, it 
is important to check how strongly the BD effect perturbs 
measured parameters of such binary systems. 



3 MEASURABILITY OF MASSES AND 

ORBITAL INCLINATION IN THE EDGE-ON 
BINARY SYSTEMS 

To investigate the measurability of the BD effect we shall 
follow the linearization procedure for evaluation of the pul- 
sar parameters, developed by Blandford & Teukolsky (1976), 
and improved by Damour & Deruelle (1986). The procedure 
is based upon the, so called, "periodic approximation" of 
differential timing formula (see below). The idea is to intro- 
duce instead of the total set of orbital and spin parameters 
of the pulsar in the timing algorithm a convenient set of 
independent parameters, which absorb all the secular time 
dependencies in the form of polynomials of time t. Coeffi- 
cients of the differential timing formula are in such case only 
periodic functions of time. 

Let the observed rotational phase of the pulsar be N b s - 
Then we can improve the initial estimates of any pulsar's 
parameter ef^ using the residuals of the pulsar's rotational 
phase R{t k ) = N b s (t k ) - N ca i(ef' ,t k ), where the numerical 
values of N ca i are computed at the moments of observed 



where the variations of the adjusted parameters in formula 
@ are: 



5e2 
6e3 

Se s 
Sep, 

5e ? 



s r 

— (SN +t- 8v v + — 



8 a 

*C9 + 7) 
-8a - t ■ 

5e 
8m c 



■ 8i> p + ■ 



8n 



(12) 



I1 1 c 

<5sinz 



with a = x sin lo and = VI -e 2 a; cos tt;; a is the orbital 
reference phase; n = 2n /Pb, Pb is the orbital period; 7 is the 
gravitational red shift plus transverse Doppler effect. Let us 
point out that to get relative accuracies in the determination 
of the mass m c and sini we have used the substitutions for 



the parameters of the type m c 

sin 



(1 + fi) and sin i 



(1 + £), where m! an d sini' ' are supposed to be 
fixed and equal to the current estimated values for these 
parameters. 

The functions Bi = dN(t)/dei are the partial deriva- 
tives of the pulsar rotational phase N(t) with respect to the 
e,:-th adjusted parameter of the pulsar, and we neglect the 
relativistic post-Keplerian parameters 8 r , 8e, k and all hid- 
den parameters (see for explanation of the physical meaning 
of the omited parameters the papers of Damour & Deruelle 
(1986) and Kopeikin (1994)): 



B 2 
B 3 

Bi 

B 5 
B 6 



cos u — e, 
sin-u, 

a sin u — P cos u 
1 — e cos u ' 

— a — B4 sin u, 
Bu, = As + As, 



B 7 = B e = L~ 1 s 



(2r + A s )(l - ecosii) 
x sin(tj + A e ). 



sC 



(13) 
(14) 
(15) 

(16) 

(17) 
(18) 

(19) 



After the first step of the least-squares fit of the residu- 
als R(t), the obtained corrections for the pulsar parameters 
Sef 1 ^ are added to their initial values = e' ' + ^e' ' and 
are used again at the second step of the iteration proce- 
dure. The process is repeated untill the limited values of the 
estimated parameters are obtained. 
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A similar model has been used by Ryba & Taylor (1991) 
and Kaspi et al. (1994) in the measurement of the param- 
eters of the PSR B1855+09. However, these authors have 
relied upon the algorithm developed by Damour & Deruelle 
(1986). The main difference between our differential timing 
formula and the one of Damour & Deruelle (1986) is in the 
expressions for the partial derivatives Bq and B-j taking into 
account the BD effect in the explicit form. It is important to 
note that omitting the BD effect from the timing algorithm 
will result, in principle, in incorrect estimates for the rela- 
tivistic post-Keplerian parameters r and s. This fact must be 
take into account for unambigious testing of the alternative 
gravity theories in the strong field regime. 

Using our improved expressions ( |l8| ) and ( p^ ) for these 
functions one can examine how the BD effect will manifest 
itself in the fitted values of the pulsar parameters /i and £. 
To investigate this question we take numerical simulations 
using as an example imaginary binary systems having the 
same parameters as the pulsars B1855+09 and B1913+16 
exept for the inclination angle i of the orbit. The inclination 
is varied in the range from 86.5° to 90°. Let us suppose 
that the standard deviation of the mock arrival times is e (a 
constant measured in fisec), and covariance matrix dj for 
the fitting method is e 2 times the inverse of the 7x7 matrix: 



Figure 2. The values of logio(cr M ) =logio(<r mc /m c ) vs. (cosi) 
for measurement of the pulsar's companion mass in the two tim- 
ing models for binary pulsars like PSR B1855+09. The dotted 
curve (S) represents the values <t m being expected from the data 
fit according to the Damour &: Deruelle (1986) timing algorithm, 
and the two solid curves represent the values cr^ being expected 
from the data fit with the gravitational bending delay taken into 
account. The assumed values for the angles n = —85°; A = 0.1°, 
and 1° correspond to the curves labeled as (0.1), and (1) respec- 
tively. The long tic on the bottom horizontal axis denotes the 
value of sin i = 0.9992 presented in (Kaspi et al. 1994). 



Bij — Bi(t k ) ■ Bj(t k ). 



(20) 



Here n is the total number of the TOA measurements and 
the functions Bi (i — 1, 2, ...7) are given by equations ( |l3| ) - 
(p"9]). Then our variance estimates for the evaluated param- 
eters ej are: 



da 



(21) 



where the errors a ei = y'bij of the parameters are computed 
by using the diagonal elements of the inverse matrix bij = 
By 1 . We also suppose in our numerical simulations that the 
observational points are equally distributed over the pulsar's 
orbital phase and their total number is equal to 1000. Two 
different timing models have been compared, namely, one of 
Damour & Deruelle (1986) without the BD effect, and that 
outlined above, which includes the BD effect. 

First of all, we took calculations for binary systems hav- 
ing orbital parameters like the PSR B1855+09. Figs 2 and 
3 illustrate the behavior of the fractional uncertainties 
= <Jm a /m c and — ti s i n i / sin i plotted as functions of cos i 
for these two algorithms respectively, 

where we fixed the value r\ = —85°, and supposed the incli- 
nation angle i ranges from 86.5° to 90°. Our assumed values 
for A = 0°.l, and 1° were chosen for illustrative purposes. 
They correspond for the case of the PSR B 1855+09 system 
to the values of the parameter D = —18ns, and —1.8ns, re- 
spectively. Note that the numerical values for the parameter 
C depend on the inclination angle, and have been computed 
from the relation C = — Dcotncosi following from equa- 
tions ( and (px|). This relation between C and D has been 
used in equations (|l^), ( |l^ ) for drawing the plots. In fact, 
the values of C for the chosen set of D and sini lie in the 
range ~0. Ins. Replacing the positive sign of the parame- 
ter D with the negative one with the same numerical values 



Figure 3. The same as Fig. 2, 
values are logio((T^) =logio(<T a i r 



except for only that the plotted 
il sini) vs. (cosi ). 



(that corresponds to the positive value of angle r\ — 85°) 
does not change the Figures 2 and 3 perceptibly. 

The same procedure has been applied to estimate the 
measurability of the "range" and "shape" parameters in the 
case of binary system PSR B1913+16. In such systems, with 
more massive pulsar companions, one can expect to measure 
the BD effect more easily. We took four numerical values 
for A = 0°.1,1°, 10°, and 80°, that corresponds to D = 
— 15721ns, —1572ns, —158ns, and —28ns. The parameter C 
lies in the range -f- 0.2ns. Results of the calculations are 
shown in the figures 4 and 5. 

As one can see from the plots (2) and (3), the inclusion 
of the effect of the gravitational bending delay is of no prac- 
tical consequence for reducing the relative errors in the mass 
and sini determinations in such systems as the binary pul- 
sars like PSR B1855+09 when cosi > 0.01 and A > 1°. As 
for the binary pulsars like PSR B1913+16, possessing more 
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Figure 4. The values of logio(<r M ) =logio(<x mc /m c ) vs - (cosi) 
for measurement of the pulsar's companion mass in the two tim- 
ing models for binary pulsars like PSR B1913+16. The dot- 
ted curve (S) represents the values er^ being expected from the 
data fit according to the Damour & Deruelle (1986) timing al- 
gorithm, and the four solid curves represent the values oy, being 
expected from the data fit with the gravitational bending delay 
taken into account. The assumed values for the angles r\ = —85°; 
A = 0°.1,1°,10°, and 80° correspond to the curves labeled as 
(0.1), (1), (10), and (80) respectively. Longitude of the perias- 
tron at the initial epoch uig = 180° . 



4 MEASURABILITY OF THE 

GRAVITATIONAL DEFLECTION OF LIGHT 
PARAMETERS IN THE PSR B1855+09 

The numerical simulations made in the previous paragraph 
show that the amplitude of the BD effect is, probably, too 
small to be detected with confidence in the PSR B1855+09 
since the real value of cos i = 0.04. Nevertheless, in the same 
spirit of "looking for zero to be sure it is there", we have 
tried to find it using the 7-yrs. observational data set for 
PSR B1855+09 obtained by J. Taylor and collaborators at 
the Arecibo observatory (Kaspi et at, 1994). To process the 
observational data we apply the timing program TIMAPR 
described in (Doroshenko & Kopeikin 1990). It has been gen- 
eralized for timing of binary pulsars in accordance with the 
timing algorithm described in Sec. 3 with a secular drift of all 
relevant parameters. The celestial reference frame was based 
upon the DE200/LE200 ephemerides of the Jet Propulsion 
Laboratory (Standish 1982) and the relativistic time scales 
transformations in the Solar system developed by Fairhead 
& Bretagnon (1989) are used in our fitting procedure (see the 
paper of Brumberg & Kopeikin (1990) for more theoretical 
details concerning the origin and structure of these trans- 
formation). The polar motion corrections were introduced 
using the "raw values of X,Y, universal time for every 5 
days" refered to the 1979 BIH system (BIH Annual Report) 
and the "normal values of the Earth orientation parameters 
at 5-day intervals" refered to the EOP (IERS) C02 system 
(IERS Annual Report). 

The general concept of measurability of pulsar parame- 
ters is described in the papers of Taylor & Weisberg (1989) 
and Damour & Taylor (1992), and we follow the outlined 
procedure considering now C and D in the formula (M) as 
the new independent fitted parameters. It changes the for- 
mula to: 



R(t) = ^ Bi ■ Se t 



(22) 



Figure 5. The same as Fig. 4, except for only that the plotted 
values are logio(cr^) = logio(cr s i ni / sini) vs. (cosi ). 



massive pulsar's companion, the BD effect has no influence 
when cosi > 0.01 and A > 80° . From an evolutionary point 
of view it is expected that for the pulsars under consider- 
ation A ~ i ~ 90° and the only real hope to measure the 
influence of the BD effect is in this case when cos i < 0.003. 
This value is close enough to that anticipated by Epstein 
(1977). It should be noted, however, that these estimates 
are true only for the total number of "observational points" 
being equal to 1000. Taking into account many more points 
or concentrating observations near the moment of superior 
conjunction makes the obtained estimates less conservative. 



where parameters ei, es are the same as in the equations 
(14), {e i; 6 < % < 9} = {r,s,C,D}. Functions B % = B % for 
i = 1,2, 3, 4, 5 are taken from the equations (|l^)-([l7|) and: 

B 6 = B r = -21nL, (23) 

B 7 = B 3 = L~ 1 (2r + A B )(1- ecosu)sm(tJ + A e ), (24) 

B$ ee Be = L 1 sin(cj + A e ), (25) 

B 9 ee B D = L' 1 cos(ui + A e ). (26) 

To determine parameters r, s, C and D we used the 
least-squares procedure of minimizing the statistic: 



•s^\N{U)-Ni 



OiV v 



(27) 



where Oi is the estimated uncertainty of the TOA t%. The 
statistic x 2 measures the goodness of fit of the timing model. 
However, direct application of the least-squares method with 
all parameters considered to be varied freely fails to deter- 
mine the parameters C and D because of their smallness. 
Moreover, the parameter C includes the additional, small 
factor cosi = 0.04, and, for this reason, has been assumed 
to be equal to zero identically. Then, to extract some valu- 
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Figure 6. The smooth curve represents values of Ax 2 obtained 
in data fit with the fixed values of the r and s relativistic post- 
Keplerian parameters obtained by Kaspi et al. (1994). The levels 
of confidence, corresponding to la, 2a and 3<r determination of 
parameter D are plotted with dotted lines. 



measurability of the magnitude r and shape s parameters 
of the Shapiro time delay as well as the BD parameters C 
and D has been examined. Processing of the observational 
data for the binary pulsar B1855+09 was rendered in an ef- 
fort to find some indication on the presence of effect of the 
gravitational deflection of light in this system. 

Taking into account the effect of the gravitational de- 
flection of light in the edge-on binary pulsars allows us to 
make more precise experimental tests of alternative theories 
of gravity in the strong gravitational field regime when the 
cosine of the orbital inclination to the plane of the sky is 
less than 0.003. This rather conservative estimate can be re- 
laxed if the rotational axis of the pulsar is close enough to 
the line of sight and/or the total number of equally spaced 
observations exceeds 1000. 

After this paper was submitted for publication, we were 
informed about a paper by Goicoechea et al. where the BD 
effect is discussed for the case of binary pulsars and X-ray 
binaries on circular orbits. 



able information, we investigate the structure of the func- 
tion Ax 2 (ei) = X 2 ( e i) ~ Xo( e i St ) in the neighborhood of the 
global minimum Xo in the multi-dimensional space of the fit- 
ted parameters (Taylor & Weisberg ,1989; Ryba and Taylor 
,1991; Damour & Taylor ,1992). 

We have calculated the functional dependence of x 2 
statistic near the region of its global minimum with the 
minimization of the x 2 by the fit of the astrometric 
{iVo, a, 5, Ha, /«}, spin {v p , v v , i> p }, orbital Keplerian {x — 
a sin i, e, To, to, Pi,}, and post-Keplerian {F%,e, u), } parame- 
ters of the pulsar with the parallax n, dispersion measure 
DM , and post-Keplerian parameters r and s held fixed at 
the unique values reported by Kaspi et al. (1994). In this 
procedure the parameter D was held fixed at a number of 
different points from -40 ns to 40 ns. The results of the fitting 
procedure are shown in Fig. 6. 

As one can see from this plot the corresponding curve 
has a parabolic shape, and the numerical value of the pa- 
rameter D = 12lg° nanosecond is not zero with the 68 % 
confidence ( lcr ) limits. The values of other astrometric, 
spin, Keplerian, and post-Keplerian parameters for the PSR 
B1855+09 are in a perfect ( lcr ) agreement with the values 
reported by Kaspi et al. (1994). Although we can not say 
anything with the adopted ( 3cr ) confidence about the real 
numerical value of the BD effect, the parabolic structure of 
the x 2 near region of its global minimum (with respect to 
the parameter D) indicates that the effect in question can 
be, perhaps, measured in the future. 



5 CONCLUSION 

We have derived the improved timing formula (^) for binary 
pulsars taking into account the relativistic effect of light de- 
flection in the gravitational field of the pulsar's companion. 
It has been shown the effect is perceptible only for the pul- 
sars with nearly edge-on orbits. The measurement of the 
effect in question permits, in principle, the determination of 
the position angle r) of the angular velocity vector of the pul- 
sar although it fails to determine direction of the rotation. 
The influence of the gravitational deflection of light on the 
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ABSTRACT 

An improved formula for the timing of binary pulsars that accounts for the relativistic 
deflection of light in the gravitational field of the pulsar's companion is presented, and 
the measurability of this effect together with its variance estimates are discussed. The 
deflection of the pulsar's beam trajectory in the gravitational field of its companion 
leads to variation in the pulsar's rotational phase. This variation appears as the narrow 
sharp growth of the magnitude of the post-fit pulsar timing residuals in the vicinity of 
the moment of the superior conjunction of the pulsar with its companion. In contrast 
to the relativistic Shapiro effect the amplitude of the effect of gravitational deflection 
of the pulsar radio beam has two peaks with opposite signs, which become sharper as 
the inclination i of the pulsar's orbit approaches 90°. The effect under consideration 
influences the estimation of parameters of the relativistic Shapiro effect in the binary 
pulsars with nearly edgewise orbits. Its inclusion in the fitting procedure provide more 
careful measurement of sine of the orbital inclination i as well as the masses of the pulsar 
m p and its companion m c . This permits an improved testing of alternative theories of 
gravity in the strong field regime. Moreover, the measurement of this effect permits 
independent geometric constraints on the position of the pulsar rotational axis in space. 
The effect of the gravitational deflection of light has been numerically investigated for 
binary pulsars with nearly edgewise orbits. It is shown that the effect is observed in 
general only when cos i < 0.003. This estimate becomes less restrictive as the pulsar's 
spin axis approaches the line of the sight. 

Key words: gravitation - pulsars - relativity - stars: neutron - stars: individual (PSRs 
B1855+09, B1913+16, B1534+12) 



1 INTRODUCTION 

High-precision timing observations of binary pulsars on rel- 
ativistic orbits provide unique opportunities for testing of 
General Relativity (GR) through two mechanisms: the emis- 
sion of gravitational radiation and the presence of strong- 
field effects. The effect of emission of gravitational waves 
have been confirmed using the PSR B1913+16 long term 
data with an accuracy better than 0.5% (Taylor & Weis- 
berg 1989; Damour &Taylor 1991). This was accomplished 
by the measurements of three relativistic effects: 1) the 
rate of apsidal motion Co, 2) the Doppler and gravitational 
red shifts 7 of pulsar rotational frequency, and 3) the sys- 
tem orbital period decay Pb, which is due to the effect of 
the gravitational waves emission, which has been rigorously 
justified from the mathematical point of view by Damour 

* On leave from Astro Space Centre of P.N.Lebedev Physical 
Institute, Leninskii prospect 53, Moscow 117924, Russia 



(1983a, 1983b) and, independently, by Grishchuk & Kopeikin 
(1983,1986), as well as Schafer (1985). However, timing of 
PSR B1913+16 is not sufficiently sensitive to the strong 
gravitational field effects. And, indeed, the theoretical inves- 
tigations of Damour and Esposito-Farese (1992) have shown 
that some alternative to the GR theories of gravity, which 
include the additional scalar fields, also fit the above men- 
tioned Co — 7 — Pb test, and, hence, can not be distinguished 
from the GR theory. Realizing this fact has led to the devel- 
opment of the parameterized post-Keplerian (PPK) formal- 
ism, which was designed to extract the maximum of possible 
information from the high-precise observational data for bi- 
nary pulsars (Damour & Taylor 1992). The application of 
this methodology to the four pulsars B1913+16, B1855+09, 
B1534+12 and B1953+29 has already ruled out a number 
of the alternative gravity theories (Taylor et al. 1992). 

The PPK formalism is especially effective in discrim- 
inating alternative theories of gravity in the strong field 
regime when there is a possibility of independently measur- 
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ing two post-Keplerian parameters chaacterizing the range 
r = Gm c /c 3 and the shape s = sin i of the famous Shapiro 
delay in propagation of radio signals of pulsar in the gravi- 
tational field of its companion. The most favorable case for 
the measurements of r and s is when the pulsar's orbit is 
inclined to the plane of the sky at the angle i close to 90°. 
Timing observations of the "non-relativistic" binary pulsar 
B1855+09 by Ryba & Taylor (1991), and by Kaspi et al. 
(1993) have provided the very first such independent mea- 
surements of the parameters r and s that allowed the deter- 
mination of the neutron star mass with a precision better 
than about 18%. It is worthwhile to stress that such mass 
determination does not require a knowledge of the rate of 
relativistic apsidal motion Co and the red shift parameter 7 
as it takes place in the case of the PSR B1913+16 system. 

At the same time when a binary pulsar orbit is viewed 
nearly edge-on (as for B1855+09 and PSR B1534+12) the 
complementary effect in the propagation of the pulsar's ra- 
dio signal may be important. This is due to the well known 
effect of the deflection of light rays in the gravitational field 
of a self-gravitating body (in the gravitational field of the 
pulsar's companion in the event of the binary pulsar system). 
Due to the gravitational deviation of the pulsar's beam from 
a straight line, the observer will see the pulsar's pulse only 
when the pulsar is rotated by an angle compensating this 
deviation. In contrast to the relativistic Shapiro effect, the 
effect of gravitational deflection of light directly contributes 
to the rotational phase of the pulsar (as, for instance, the 
effect of aberration of pulsar's beam does), but this is so 
small that it can be misinterpreted as an additional delay 
in the time-of-arrivals (TOA) of the pulsar pulses. Similar 
to the Shapiro delay, the deflection of the pulsar's beam by 
the gravitational field of its companion manifests itself in 
the timing as the rapid, sharp growth of the magnitude of 
the post-fit residuals of TOA on a short time interval in the 
vicinity to the moment of superior conjunction of the pulsar 
and its companion. The effect in question is superimposed 
on the Shapiro effect, and, in principle, should be explicitly 
taken into account to avoid an incorrect determination of 
the parameters r and s. 

Another important feature of the effect of gravitational 
deflection of light in a binary pulsar is a strong dependence 
of the shape of TOA residuals on the spatial orientation of 
the pulsar's rotation axis (see the eqs. (7), (9), (10) in the 
text below). Thus, the measurement of this effect would help 
to get further constrains on the position of the pulsar's spin 
independently of the polarization observations. Let us em- 
phasize that the effect of the relativistic aberration caused 
by the orbital motion of pulsar in binary systems (Smarr 
& Blandford 1976) also strongly depends on the spin's ori- 
entation, but can not be observed directly due to the same 
time-dependence on the pulsar's orbital phase as the classi- 
cal Roemer effect (Damour & Deruelle 1986). 

Epstein (1977) was, probably, the first who called at- 
tention to the phenomenon of deflection of the pulsar ra- 
dio beam in the gravitational field of its companion when 
cos i < 10 -3 . However, he did not make any calculations 
confirming this estimate. Narayan et al. (1991) also men- 
tioned this effect as a useful tool for investigation of black 
hole - neutron star binary systems and verification of the 
existence of the black hole to an unparalleled degree of cer- 
tainty. Schneider (1989, 1990) calculated the corresponding 



gravitational deflection time delay and lensing amplification 
factor of the pulse intensity in the edge-on binary pulsar, 
but, unfortunately, his presentation is not easily followed. 

In the present paper we try to solve the problem more 
simply, giving both a complete theoretical treatment and a 
discussion of the measurability of the effect of the gravita- 
tional deflection of radio beam in binary pulsars which we 
call briefly the bending delay (BD) effect. It is shown how 
to include this effect in the fitting procedure of TOA for the 
evaluation of the pulsar parameters of binary pulsars whose 
orbital inclinations to the plane of the sky are extremely 
close to 90° . 

In Sec. 2 of the paper we present an improved formula 
for the timing of binary pulsars with the BD effect taken into 
account. In Sec. 3 we investigate, using numerical simula- 
tions, the measurability of the masses and orbital inclination 
in binary systems when the BD effect is taken into account. 
A mathematical procedure for the direct measurement of 
parameters of the BD effect is outlined in Sec. 4 and applied 
to process the observational data for PSR B1855+09. In Sec. 
5 we present a summary of our results and conclusions. 



2 TIMING FORMULA 

The self-consistent mathematical derivation of a timing for- 
mula for the calculation of pulsar rotational phase should 
be based upon the relativistic theory of astronomical refer- 
ence systems together with the matching procedure devel- 
oped previously by Kopeikin (1988), Brumberg & Kopeikin 
(1989a, b) and Brumberg (1991). The exact timing formula 
must include the kinematic effects of the galactic differen- 
tial rotation, orbital and precessional motions of the pulsar, 
relativistic effect of aberration, quadratic Doppler and red 
shift effects, the Shapiro delay, as well as the deflection of 
the pulsar beam in the gravitational field of its compan- 
ion. In the present paper we ignore the effects of differential 
galactic acceleration and precessional motion of the pulsar 
rotational axis suggesting for simplicity that vectors of the 
angular velocity Q and spin S of the pulsar are coincident. 
The timing formula obtained is a slightly modified version 
of that elaborated by Damour & Deruelle (1986). Some ad- 
ditional theoretical aspects of the derivation of the timing 
formula are discussed by Klioner & Kopeikin (1994) and 
Kopeikin (1994, 1995). 

Let ri and A be the polar coordinates of the pulsar's 
spin axis 5, where 77 is the longitude of the pulsar's spin 
axis in the plane of the sky measured from the ascending 
node of the pulsar's orbital plane (0° < 77 < 360°), and A 
is the angle between the pulsar's spin axis and the line of 
sight pointing from the terrestrial observer toward the pulsar 
(0° < A < 180°). The time evolution of the pulsar rotational 
phase 4>{T) is given by the equation: 

0(T)/2tt = N{T) = N + v P T + h p T 2 + h p T 3 , (1) 

2 

where No is the initial phase; v p , v p , v p denote the pulsar spin 
frequency and its time derivatives shifted by the Doppler 
factor, and T is the pulsar proper time related to the solar- 
system barycentric arrival time t according to the relation- 
ship: 

t-t = T + A R (T) + A E (T) + As(T) + A A (T) + A B (T).(2) 



Relativistic Effect of Gravitational Deflection of Light in Binary Pulsars 3 



Here io is the initial barycentric epoch of observation; Ah, 
Ae and A^i correspond to the well-known Roemer, Ein- 
stein and aberration time delays, whose exact mathematical 
expressions can be found in the paper of Damour & Tay- 
lor (1992); As is the Shapiro delay of radio pulses in the 
gravitational field of companion of the pulsar (Epstein 1977; 
Damour & Taylor 1992): 

A s = -2r In L, (3) 

L = 1 — e cos u — s{sin cj(cos u — e) + (l — e 2 ) 1 '' 2 cos u sin u}, (4) 

which depends on the just introduced post-Keplerian pa- 
rameters r and s as well as the classical orbital parameters: 
eccentric anomaly u, eccentricity e, and longitude of the pe- 
riastron ui. 

Ab in the equation (2) is the additional bending delay 
caused by the deflection of the pulsar beam in the gravi- 
tational field of its companion. It has been calculated from 
the condition of emission of the pulsar's pulse toward the 
observer (Kopeikin 1992). This condition is expressed in the 
pulsar's proper reference frame by the equation: 



(JV ■ (fi A b)) = 0, 



(5) 



which defines a relation between the number of pulse, which 
has been emitted, and the moment of the emission reckoned 
in the pulsar proper time scale T. Here N denotes the spatial 
components of isotropic unit vector tangent to the trajectory 
of the emitted pulse, and directed toward the terrestrial ob- 
server; b is the unit vector in the fiduciary direction of maxi- 
mal pulsar radio emission (it is not necessarily the direction 
of magnetic dipole); Q is the vector of the instantaneous an- 
gular velocity of the pulsar; the signs "■" and "A" denote 
the Euclidean scalar and cross products. We have assumed 
implicitly that the model of the conical rotating beacon for 
the description of the pulsar's emission (Lyne & Manchester 
1988) is true. However, we would like to point out that cal- 
culations being described in this paragraph are still valid in 
the case of a deviation of the polar beam geometry from the 
conical one as it is, for instance, in the model of Narayan & 
Vivekanand (1983). 

Vector N does not keep a constant spatial orientation 
from the point of view of an external observer at rest but 
has small (caused by both classical and relativistic effects) 
spatial variations with respect to the constant unit vector n 
referred to the rest frame of the binary pulsar barycenter. 
To show this we use the simplified version of the relativistic 
relationship between the spatial components of the vectors 
N and n derived by Klioner & Kopeikin (1992) in the case 
of the relativistic N-body problem: 



(6) 



N = ft — ^(™ A (f p A ft)) + -(ft A (v p A ft)) 
2Gm c (ft A (tr A ft)) 

C 2 TR TR + (ft ■ tr) 

to replace in the equation (5) the time dependent compo- 
nents N by the ones of the vector n. The complete relation 
between N and n (Klioner & Kopeikin 1992) can be easily 
used for investigations of the second order aberration effects 
and physical parameters of a rotating black hole being con- 
sidered as a possible companion of a pulsar (Narayan et al. 
1991). 

In formula (6), c is the speed of light, r p and tr are 



the radius- vectors pointing respectively from the barycenter 
of the binary system and the center-of-mass of the pulsar's 
companion to the pulsar (do not confuse these vectors with 
the notation r for the parameter characterizing amplitude 
of the Shapiro delay), v p = dfp/dT is the orbital velocity of 
the pulsar, and R is the distance between the binary pul- 
sar and the solar system barycenters. Let us note that it is 
possible to equate the vector —n to the unit vector K (that 
is n = —K) pointing from the Earth toward the pulsar sys- 
tem as it is usually done, though such procedure requires an 
appropriate justification from a theoretical point of view. 

After substituting expression (6) into equation (5), ex- 
pressing vectors N , Q, and n through the Euler variables, 
and solving the obtained trigonometric equation with re- 
spect to the fast Euler variable as described in (Doroshenko 
et al. 1995 we get the counting formula (1) with additional 
phase corrections due to the classical and relativistic terms 
in the relationship (6). To be specific, the second term in the 
right-hand side of this formula gives an additional but neg- 
ligibly small classical parallax delay in the rotational phase 
of the pulsar (see, however, the paper of Kopeikin (1995) 
where the influence of the orbital parallax effect on the 
propagation time of pulsar radio signals is shown to be de- 
tectable); the third term in the right-hand side of (6) leads 
to the aberration delay A^(T) depending on two directly 
unobservable parameters A and B introduced by Damour & 
Deruelle (1986), and the very last term defines after straight- 
forward calculations the gravitational bending delay being 
under consideration: 



A B = L 1 (C sin(w + A e ) + D cos(w + A e )) . 



(7) 



Here the true orbital anomaly A e is related to u by the well- 
known equation: 



A e = 2 arctan 



1 + eN 1 / 2 u 

tan — 

1 - e) 2 



(8) 



clr designates the semi-major axis of the relative orbit of 
the pulsar with respect to its companion, and the constants 
C and D have been introduced as the new post-Keplerian 
parameters: 



C : 



D 



r(l - s 2 ) 1/2 cos 7] 



irVpdR 



sin A 



sin rj 



itVpCLR sin A 



(9) 



(10) 



where uh = a p (m p + m c )/m c , and a p is the major semi- 
axis of the pulsar orbit with respect to the barycenter of 
the binary system. It is interesting to note that the final 
expression for the BD effect, formula (7), does not include 
any dependence on the spatial orientation of vector b. The 
same is true for the aberration delay A^. 

For a better understanding of influence of the BD effect 
on the TOA of the pulsar pulses let us compare the behavior 
of the functions As and Ab by inserting parameters of the 
system PSR B1855+09 into the expressions (3) and (7). The 
explicit dependence on time of the Shapiro (eq. (3)) and 
BD (eq. (7)) effects are shown in Fig.l , which has been 
drawn using observational parameters of the PSR B1855+09 
reported by Kaspi et al. (1994). 

For illustrative purposes we have chosen the artificially 
small numerical value for A = 0°.005, and two different val- 
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Figure 1. The illustrative example of the TOA residuals for PSR 
B1855+09 caused by the Shapiro delay and deflection of the pul- 
sar beam in the gravitational field of its companion - smooth or 
dotted curves. For more explanations see the text. 



ues for the angle: i) = — 85° (solid curve, D > 0, C > 0), 
and ri = 85° (dotted curve, D < 0, C > 0) to depict the 
BD effect in the most clear form. One can see from Fig.l 
that the modulation of TOA due to the BD effect is com- 
pletely different and quite sharper than that caused by the 
Shapiro delay. The total observational timing effect will be 
the sum of the Shapiro and BD delay curves. It is important 
to note that the asymmetry in the curve describing the BD 
effect does not depend on the value of the angle A (since 
<sin A < 1, when 0° < A < 180°) and is sensitive to the 
numerical value of the angle i) only. Therefore, we can say 
nothing from the investigation of the BD effect whether the 
pulsar spin axis S is along the direction to observer or in 
the opposite one. 

The considered example shows that there is in principal 
a possibility to detect the BD effect in the binary pulsars 
with nearly edgewise orbits. However, this effect was never 
included in existing pulsar timing algorithms. Therefore, it 
is important to check how strongly the BD effect perturbs 
measured parameters of such binary systems. 



3 MEASURABILITY OF MASSES AND 

ORBITAL INCLINATION IN THE EDGE-ON 
BINARY SYSTEMS 

To investigate the measurability of the BD effect we shall 
follow the linearization procedure for evaluation of the pul- 
sar parameters, developed by Blandford & Teukolsky (1976), 
and improved by Damour & Deruelle (1986). The procedure 
is based upon the, so called, "periodic approximation" of 
differential timing formula (see below). The idea is to intro- 
duce instead of the total set of orbital and spin parameters 
of the pulsar in the timing algorithm a convenient set of 
independent parameters, which absorb all the secular time 
dependencies in the form of polynomials of time t. Coeffi- 
cients of the differential timing formula are in such case only 
periodic functions of time. 

Let the observed rotational phase of the pulsar be Nobs- 
Then we can improve the initial estimates of any pulsar's 
parameter e'°' using the residuals of the pulsar's rotational 
phase R(t k ) = Nob^-Ncai^Jk), where the numerical 
values of N ca i are computed at the moments of observed 



arrival times tk with the help of the equations (1) - (4) and 
(7) - (10). Assuming that the differences Sei = ei — e'°' 
between the (yet unknown) true and initial numerical values 
of the parameters are small, one can represent an expression 
for the total derivative of the phase residuals (the differential 
timing formula) in the form: 



R(t) = ^ B i ■ Se i 



(11) 



where the variations of the adjusted parameters in formula 



(11) are: 




8ei = 


- — (6N + t 
i/ p 


6e 2 = 


6a 


6e 3 = 


ffG0+7) 


8ei = 


— 6cr — t ■ 6n 


6es = 


6e 




6m z 


de^ = 






m z 




6 sin i 


8e( = 


sin i ' 



(12) 



with a = x sin w and /3 = \/l — e 2 xcoslj; a is the orbital 
reference phase; n = 2ir/Pb, Pb is the orbital period; 7 is the 
gravitational red shift plus transverse Doppler effect. Let us 
point out that to get relative accuracies in the determination 
of the mass m c and sini we have used the substitutions for 
the parameters of the type m c = m' ' ■ (1 + /1) and sin i = 
sini' ' ■ (1 + £), where m' ' and sini' ' are supposed to be 
fixed and equal to the current estimated values for these 
parameters. 

The functions Bi = 8N(t)/ dei are the partial deriva- 
tives of the pulsar rotational phase N(t) with respect to the 
ei-th adjusted parameter of the pulsar, and we neglect the 
relativistic post-Keplerian parameters 6 T , 6e, k and all hid- 
den parameters (see for explanation of the physical meaning 
of the omited parameters the papers of Damour & Deruelle 
(1986) and Kopeikin (1994)): 



Si = 1, 

B2 = cos u — e, 
B3 = sin u, 

a sin u — f3 cos u 

-D4 = " , 

1 — e cos 11 
Be, = —a — Bi sin u, 

B 6 = Bu. = As + As, 



(13) 
(14) 
(15) 

(16) 

(17) 
(18) 



B 1 = B i = |(2t- + A B )(1 - ecosa) - — -j 

xsin(w + J 4 e ). (19) 

After the first step of the least-squares fit of the residu- 
als R(t), the obtained corrections for the pulsar parameters 
(5e'°' are added to their initial values e' 1 ' = e'°' + 5e'°' and 
are used again at the second step of the iteration proce- 
dure. The process is repeated untill the limited values of the 
estimated parameters are obtained. 



Relativistic Effect of Gravitational Deflection of Light in Binary Pulsars 5 



A similar model has been used by Ryba & Taylor (1991) 
and Kaspi et al. (1994) in the measurement of the param- 
eters of the PSR B1855+09. However, these authors have 
relied upon the algorithm developed by Damour & Deruelle 
(1986). The main difference between our differential timing 
formula and the one of Damour & Deruelle (1986) is in the 
expressions for the partial derivatives B$ and Br taking into 
account the BD effect in the explicit form. It is important to 
note that omitting the BD effect from the timing algorithm 
will result, in principle, in incorrect estimates for the rela- 
tivistic post-Keplerian parameters r and s. This fact must be 
take into account for unambigious testing of the alternative 
gravity theories in the strong field regime. 

Using our improved expressions (18) and (19) for these 
functions one can examine how the BD effect will manifest 
itself in the fitted values of the pulsar parameters fi and £. 
To investigate this question we take numerical simulations 
using as an example imaginary binary systems having the 
same parameters as the pulsars B1855+09 and B1913+16 
exept for the inclination angle i of the orbit. The inclination 
is varied in the range from 86.5° to 90°. Let us suppose 
that the standard deviation of the mock arrival times is e (a 
constant measured in /zsec), and covariance matrix Cij for 
the fitting method is e 2 times the inverse of the 7x7 matrix: 

n 

Bij =^Bi(i fc ). Bj(t k ). (20) 
fc=i 

Here n is the total number of the TOA measurements and 
the functions Bi (i = 1, 2, ...7) are given by equations (13) - 
(19). Then our variance estimates for the evaluated param- 
eters ei are: 

Sei = a ei ■ , (21) 

V n 

where the errors a ei = ^Jhij of the parameters are computed 
by using the diagonal elements of the inverse matrix bij = 
Bf. 1 . We also suppose in our numerical simulations that the 
observational points are equally distributed over the pulsar's 
orbital phase and their total number is equal to 1000. Two 
different timing models have been compared, namely, one of 
Damour & Deruelle (1986) without the BD effect, and that 
outlined above, which includes the BD effect. 

First of all, we took calculations for binary systems hav- 
ing orbital parameters like the PSR B1855+09. Figs 2 and 
3 illustrate the behavior of the fractional uncertainties cr^ 
= <7 mc I m c and = ir s i n i / sin i plotted as functions of cos i 
for these two algorithms respectively, 

where we fixed the value n = —85°, and supposed the incli- 
nation angle i ranges from 86.5° to 90°. Our assumed values 
for A = 0°.l, and 1° were chosen for illustrative purposes. 
They correspond for the case of the PSR B1855+09 system 
to the values of the parameter D = — 18ns, and — 1.8ns, re- 
spectively. Note that the numerical values for the parameter 
C depend on the inclination angle, and have been computed 
from the relation C = — i? cot t; cos i following from equa- 
tions ( 9) and (10). This relation between C and D has been 
used in equations (18), (19) for drawing the plots. In fact, 
the values of C for the chosen set of D and sin i lie in the 
range -i-O.lns. Replacing the positive sign of the parame- 
ter D with the negative one with the same numerical values 




0.01 0.02 0.03 0.04 0.05 0.06 

cos i 



Figure 2. The values of logio(o>) =logio(o" mc /m c ) vs. (cosi) for 
measurement of the pulsar's companion mass in the two timing 
models for binary pulsars like PSR B1855+09. The dotted curve 
(S) represents the values being expected from the data fit ac- 
cording to the Damour & Deruelle (1986) timing algorithm, and 
the two solid curves represent the values being expected from 
the data fit with the gravitational bending delay taken into ac- 
count. The assumed values for the angles 77 = —85°; A = 0.1°, and 
1° correspond to the curves labeled as (0.1), and (l) respectively. 
The long tic on the bottom horizontal axis denotes the value of 
sini = 0.9992 presented in (Kaspi et al. 1994). 




0.01 0.02 0.03 0.04 0.05 

cos i 



Figure 3. The same as Fig. 2, except for only that the plotted 
values are logio(<7^) =logio(<7 s i n if sini) vs. (cosi ). 

(that corresponds to the positive value of angle n = 85°) 
does not change the Figures 2 and 3 perceptibly. 

The same procedure has been applied to estimate the 
measurability of the "range" and "shape" parameters in the 
case of binary system PSR B1913+16. In such systems, with 
more massive pulsar companions, one can expect to measure 
the BD effect more easily. We took four numerical values 
for A = 0°.1,1°, 10°, and 80°, that corresponds to D = 
— 15721ns, —1572ns, — 158ns, and — 28ns. The parameter C 
lies in the range -j- 0.2ns. Results of the calculations are 
shown in the figures 4 and 5. 

As one can see from the plots (2) and (3), the inclusion 
of the effect of the gravitational bending delay is of no prac- 
tical consequence for reducing the relative errors in the mass 
and sin i determinations in such systems as the binary pul- 
sars like PSR B1855+09 when cosi > 0.01 and A > 1°. As 
for the binary pulsars like PSR B1913+16, possessing more 
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Figure 4. The values of logio(<7^) =logio(<7m c /m c ) vs. (cosi) 
for measurement of the pulsar's companion mass in the two 
timing models for binary pulsars like PSR B1913+16. The dot- 
ted curve (S) represents the values <7^ being expected from the 
data fit according to the Damour & Deruelle (1986) timing al- 
gorithm, and the four solid curves represent the values <7^ being 
expected from the data fit with the gravitational bending delay 
taken into account. The assumed values for the angles 77 = —85°; 
A = 0°.1,1 ,10°, and 80° correspond to the curves labeled as 
(0.1), (1), (10), and (80) respectively. Longitude of the perias- 
tron at the initial epoch ojq = 180°. 
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Figure 5. The same as Fig. 4, except for only that the plotted 
values are logio(<7^) = £ogio(<7 s in {/ sini) vs. (cosi ). 



massive pulsar's companion, the BD effect has no influence 
when cos i > 0.01 and A > 80° . From an evolutionary point 
of view it is expected that for the pulsars under consider- 
ation A ~ i ~ 90° and the only real hope to measure the 
influence of the BD effect is in this case when cos i < 0.003. 
This value is close enough to that anticipated by Epstein 
(1977). It should be noted, however, that these estimates 
are true only for the total number of "observational points" 
being equal to 1000. Taking into account many more points 
or concentrating observations near the moment of superior 
conjunction makes the obtained estimates less conservative. 



4 MEASURABILITY OF THE 

GRAVITATIONAL DEFLECTION OF LIGHT 
PARAMETERS IN THE PSR B1855+09 

The numerical simulations made in the previous paragraph 
show that the amplitude of the BD effect is, probably, too 
small to be detected with confidence in the PSR B1855+09 
since the real value of cos i = 0.04. Nevertheless, in the same 
spirit of "looking for zero to be sure it is there", we have 
tried to find it using the 7-yrs. observational data set for 
PSR B1855+09 obtained by J. Taylor and collaborators at 
the Arecibo observatory (Kaspi et al., 1994). To process the 
observational data we apply the timing program TIMAPR 
described in (Doroshenko & Kopeikin 1990). It has been gen- 
eralized for timing of binary pulsars in accordance with the 
timing algorithm described in Sec. 3 with a secular drift of all 
relevant parameters. The celestial reference frame was based 
upon the DE200/LE200 ephemerides of the Jet Propulsion 
Laboratory (Standish 1982) and the relativistic time scales 
transformations in the Solar system developed by Fairhead 
& Bretagnon (1989) are used in our fitting procedure (see the 
paper of Brumberg & Kopeikin (1990) for more theoretical 
details concerning the origin and structure of these trans- 
formation). The polar motion corrections were introduced 
using the "raw values of X,Y, universal time for every 5 
days" refered to the 1979 BIH system (BIH Annual Report) 
and the "normal values of the Earth orientation parameters 
at 5-day intervals" refered to the EOP (IERS) C02 system 
(IERS Annual Report). 

The general concept of measurability of pulsar parame- 
ters is described in the papers of Taylor & Weisberg (1989) 
and Damour & Taylor (1992), and we follow the outlined 
procedure considering now C and D in the formula (7) as 
the new independent fitted parameters. It changes the for- 
mula (11) to: 



R(t) = ^ Bi ■ Se { 



(22) 



where parameters ei, ...,e<g, are the same as in the equations 
(14), {e;;6 < i < 9} = {r, s,C, D}. Functions B { = B { for 
i = l,2, 3, 4, 5 are taken from the equations (13)-(17) and: 

B 6 = B r = -21ni, (23) 

B 7 = B s = L~ 1 (2r + A B )(1 - ecosu)sin(w + A e ), (24) 

B s = B c = i" 1 sin(w + A e ), (25) 

B 9 = B D = i" 1 cos(w + A e ). (26) 

To determine parameters r, s,C and D we used the 
least-squares procedure of minimizing the statistic: 



X 



E 



(27) 



where Ui is the estimated uncertainty of the TOA t{. The 
statistic % 2 measures the goodness of fit of the timing model. 
However, direct application of the least-squares method with 
all parameters considered to be varied freely fails to deter- 
mine the parameters C and D because of their smallness. 
Moreover, the parameter C includes the additional, small 
factor cos i = 0.04, and, for this reason, has been assumed 
to be equal to zero identically. Then, to extract some valu- 
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Figure 6. The smooth curve represents values of Ax 2 obtained 
in data fit with the fixed values of the r and s relativistic post- 
Keplerian parameters obtained by Kaspi et al. (1994). The levels 
of confidence, corresponding to la, 2a and 3a determination of 
parameter D are plotted with dotted lines. 



measurability of the magnitude r and shape s parameters 
of the Shapiro time delay as well as the BD parameters C 
and D has been examined. Processing of the observational 
data for the binary pulsar B1855+09 was rendered in an ef- 
fort to find some indication on the presence of effect of the 
gravitational deflection of light in this system. 

Taking into account the effect of the gravitational de- 
flection of light in the edge-on binary pulsars allows us to 
make more precise experimental tests of alternative theories 
of gravity in the strong gravitational field regime when the 
cosine of the orbital inclination to the plane of the sky is 
less than 0.003. This rather conservative estimate can be re- 
laxed if the rotational axis of the pulsar is close enough to 
the line of sight and/or the total number of equally spaced 
observations exceeds 1000. 

After this paper was submitted for publication, we were 
informed about a paper by Goicoechea et al. where the BD 
effect is discussed for the case of binary pulsars and X-ray 
binaries on circular orbits. 



able information, we investigate the structure of the func- 
tion A% 2 (ei) = % 2 ( e i) — Xo( e f s ') m the neighborhood of the 
global minimum % 2 m the multi-dimensional space of the fit- 
ted parameters (Taylor & Weisberg ,1989; Ryba and Taylor 
,1991; Damour & Taylor ,1992). 

We have calculated the functional dependence of % 2 
statistic near the region of its global minimum with the 
minimization of the % 2 by the fit of the astrometric 
{No, a, 6, Ha, Us}, spin {v p , v p , Vp], orbital Keplerian {x = 
a sin i, e, To , w, Pb], and post-Keplerian {Pj,,e,cj,} parame- 
ters of the pulsar with the parallax 7r, dispersion measure 
DM , and post-Keplerian parameters r and s held fixed at 
the unique values reported by Kaspi et al. (1994). In this 
procedure the parameter D was held fixed at a number of 
different points from -40 ns to 40 ns. The results of the fitting 
procedure are shown in Fig. 6. 

As one can see from this plot the corresponding curve 
has a parabolic shape, and the numerical value of the pa- 
rameter D = 12^g° nanosecond is not zero with the 68 % 
confidence ( la ) limits. The values of other astrometric, 
spin, Keplerian, and post-Keplerian parameters for the PSR 
B1855+09 are in a perfect ( If ) agreement with the values 
reported by Kaspi et al. (1994). Although we can not say 
anything with the adopted ( 3it ) confidence about the real 
numerical value of the BD effect, the parabolic structure of 
the % 2 near region of its global minimum (with respect to 
the parameter D) indicates that the effect in question can 
be, perhaps, measured in the future. 



5 CONCLUSION 

We have derived the improved timing formula (2) for binary 
pulsars taking into account the relativistic effect of light de- 
flection in the gravitational field of the pulsar's companion. 
It has been shown the effect is perceptible only for the pul- 
sars with nearly edge-on orbits. The measurement of the 
effect in question permits, in principle, the determination of 
the position angle i) of the angular velocity vector of the pul- 
sar although it fails to determine direction of the rotation. 
The influence of the gravitational deflection of light on the 
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